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ABSTRACT 

The  nonlinear  development  of  disturbances  in  stratified  shear  flows  (having  a  local  Richard¬ 
son  number  of  value  less  than  one  quarter)  is  considered.  Such  modes  are  initially  fast  grow¬ 
ing  but,  like  related  studies,  we  assume  that  the  viscous,  non-parallel  spreading  of  the  shear 
layer  results  in  them  evolving  in  a  linear  fashion  until  they  reach  a  position  vhere  their 
amplitudes  are  large  enough  and  their  growth  rates  have  diminished  sufficiently  so  that  am¬ 
plitude  equations  can  be  derived  using  weakly  nonlinear  and  non-equilibrium  critical-layer 
theories.  Four  different  basic  integro-differential  amplitude  equations  are  possible,  including 
one  due  to  a  novel  mechanism;  the  relevant  choice  of  amplitude  equation,  at  a  particular 
instance,  being  dependent  on  the  relative  sizes  of  the  disturbance  amplitude,  the  growth  rate 
of  the  disturbance,  its  wavenumber  and  the  viscosity  of  the  fluid.  This  richness  of  choice  of 
possible  nonlinearities  arises  mathematically  from  the  indicial  Frobenius  roots  of  the  govern¬ 
ing  linear  inviscid  equation  (the  Taylor-Goldsteiu  equation)  not,  in  general,  differing  by  an 
integer.  The  initial  nonlinear  evolution  of  a  mode  will  be  governed  by  an  integro-differential 
amplitude  equations  with  a  cubic  nonlinearity  but  the  resulting  significant  increase  in  the 
size  of  the  disturbance’s  amplitude  leads  on  to  the  next  stage  of  the  evolution  process  where 
the  evolution  of  the  mode  is  governed  by  an  integro-differential  amplitude  equations  with  a 
quintic  nonlinearity.  Continued  growth  of  the  disturbance  amplitude  is  expected  during  this 
stage,  resulting  in  the  effects  of  nonlinearity  spreading  to  outside  the  critical  level,  by  which 
time  the  flow  has  become  fully  nonlinear. 
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tract  No.  NAS1-19480  while  the  three  authors  were  in  residence  at  the  Institute  for  Computer  Applications 
in  Science  and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23681. 
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1  Introduction 


This  paper  extends  the  theories  of  Churilov  &  Shukhman  (1988)  and  Goldstein  &  Leib 
(1989)  to  the  case  of  (non-marginally)  unstable  stratified  flows.  Such  flows,  in  which  the 
vertical  density  variation  is  insufficient  to  overcome  the  strongly  destabilising  effects  of  the 
vertical  velocity  variation,  occur  naturally  in  many  situations.  The  mathematical  formulation 
of  this  problem  has  great  similarity  with,  amongst  others,  that  for  the  related  problem  of 
buoyancy-driven  instabilities  above  a  heated  plate  and  with  the,  at  first  sight,  less  physically 
related  problem  of  vortex  disturbances  driven  by  wall  curvature  occurring  in  3D  boundary 
layers.  Thus,  whilst  here  we  concentrate  solely  on  a  particular  model  stratified  shear  flow,  the 
mathematical  theory  has  much  wider  applications  to  both  geophysical  and  aerodynamical  flow 
situations. 

The  stratification  of  a  shear  flow  is  usually  characterised  by  a  physical  parameter  called  the 
Richardson  number  (herein  after  denoted  J);  it  is  a  measure  of  the  ratio  of  vertical  density 
variation  to  the  scale  of  vertical  velocity  shear.  The  Richardson  number  can  be  both  positive 
and  negative;  the  former  case  corresponding  to  lighter  fluid  lying  above  heavier  fluid,  and  the 
latter  case  corresponding  to  the  lighter  fluid  lying  below.  The  stability  of  a  stratified  shear  flow 
was  first  considered  independently  by  Taylor  (1931)  and  Goldstein  (1931).  They  concluded  that 
a  multi-layer  system  of  homogeneous  fluids  can  not  be  used  to  approximate  the  stability  of  a 
heterogeneous  fluid,  and  that  for  the  flow  models  considered,  the  flow  is  stable  for  Richardson 
numbers  greater  than  one  quarter. 

The  flow  model  adopted  in  this  paper  was  introduced  by  Drazin  (1958);  the  vertical  ve¬ 
locity  and  density  distribution  are  represented  by  hyperbolic  tangent  and  exponential  profiles, 
respectively.  This  model  enables  an  exact  solution  of  the  governing  linear  stability  equation, 
the  Taylor-Goldstein  equation,  for  neutral  disturbances  and  the  same  upper  bound  on  the 
Richardson  number  is  found.  Soon  after  Drazin’s  paper  appeared,  Miles  (1961)  presented 
several  theorems  concerning  properties  of  solutions  to  the  Taylor-Goldstein  equation.  In  par¬ 
ticular,  he  proved  the  two  important  results  that  in  general  flows  are  stable  for  J  >  1/4  and 
that  neutral  modes  are  proportional  to  just  one  of  the  associated  Frobenius  series  near  the 
critical  level.  The  continuation  of  the  linear  eigenfunctions  in  the  neighbourhood  of  the  critical 
levels  was  considered  by  Koppel  (1964). 

Several  linear  results  for  the  alternative  Holmboe  model  were  computed  by  Hazel  (1972), 
for  a  variety  of  boundary  conditions;  however,  our  interest  herein  is  in  the  nonlinear  critical- 
layer  theory  and  not  with  a  discussion  about  which  model  is  superior.  Around  the  same  time, 
nonlinear  studies  started  to  appear  eg.  Maslowe  (1972)  using  the  equilibrium  critical  layer 
theory  of  Benney  k  Bergeron  (1969).  Subsequent  nonlinear  studies  were  also  mainly  for  the 
marginal  instability  case  (0  <  1/4  —  J  1)  eg.  by  Maslowe  (1977)  who  considered  a  viscous 
critical  layer,  and  by  Brown  k  Stewartson  (1978)  who  considered  a  time-dependent  critical 
layer  for  small  time.  These  two  theories  were  unified  by  Brown,  Rosen  k  Maslowe  (1981).  The 
latter  paper  was  corrected  to  some  extent  by  Churilov  k  Shukhman  (1987);  soon  after,  Churilov 
k  Shukhman  (1988)  essentially  extended  and  corrected  the  work  of  Brown  k  Stewartson  (1978) 
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and  Churilov  &  Shukhman  (1987).  Recently,  Troitskaya  (1991)  has  considered  nonlinear  effects 
for  J  >  1/4. 

There  have  been  many  more  studies  of  (non-stratified)  shear  layers,  in  which  there  is  no 
vertical  density  variation.  The  linear  stability  of  such  flows  is,  in  general,  governed  by  the 
familiar  Rayleigh  equation,  to  which  the  Taylor-Goldstein  equation  reduces  for  zero  Richard¬ 
son  number.  The  nonlinear  stability  of  such  flows  has  received  much  attention  over  the  past 
couple  of  decades.  Benney  Sz  Bergeron  (1969)  developed  the  so-called  equilibrium  critical  layer 
theory;  here  the  mode  is  treated  as  ‘quasi-steady’  inside  the  critical  layer  as  well  as  outside 
it.  Nonlinearity  affects  the  jump  imposed  across  the  critical  layer  and  hence  leads  to  modified 
results  for  the  neutral  (equilibrium)  modes.  Haberman  (1972)  extended  the  theory  to  include 
critical  layers  where  viscosity  is  also  significant.  The  first  studies  of  non-equilibrium  critical 
layers  include  the  papers  by  Stewartson  (1978),  Warn  &  Warn  (1978),  Hickernell  (1984)  and 
Gajjar  &  Smith  (1985).  The  key  paper  by  Hickernell  (1984)  concerned  a  shear  layer  affected 
by  Coriolis  (rotational)  effects;  here  the  weakly  nonlinear  theory  leads  to  an  integro-differential 
equation  rather  than*  the  (previously)  more  familiar  Landau  equation  with  ‘polynomial’  non¬ 
linear  terms.  In  fact,  such  integro-differential  equations  result  naturally  from  non-equilibrium 
nonlinear  critical  layer  theories  when  the  shear  layer  is  coupled  with  other  physical  factors  such 
as,  for  instance,  Coriolis  effects  (eg.  Hickernell,  1984;  Shukhman,  1991);  compressibility  effects 
(eg.  Goldstein  &  Leib,  1989);  three-dimensionality  effects  (eg.  Goldstein  &  Choi,  1989;  Wu 
et  al,  1993);  and  buoyancy  effects  (eg.  Churilov  &  Shukhman,  1988).  Moreover,  in  a  related 
paper  (Blackaby,  Dando  &  Hall,  1993),  we  demonstrate  that  such  an  integro-differential  type 
of  equation  can  be  derived  to  describe  the  nonlinear  evolution  of  the  inviscid  Gortler  modes 
studied  by  Bassom  &  Hall  (1991).  However,  the  case  of  a  ‘simple’  shear  layer,  not  affected 
by  any  additional  physical  factors,  is  a  special  case  in  the  sense  that  it  does  not  lead  to  an 
integro-differential  equation;  instead,  Goldstein  &;  Leib  (1988),  found  that  the  nonlinear  evolu¬ 
tion  of  a  disturbance  was  governed  by  the  full  unsteady  nonlinear  critical-layer  equations.  This 
difference  is  due  to  the  additional  physical  factors,  of  the  former  cases,  resulting  in  stronger 
singularities  of  the  inviscid  disturbance  quantities  at  the  critical  level. 

At  first  sight,  it  appears  that  weakly  nonlinear  theories  can  only  be  usefully  applied  to 
marginally  unstable  flows;  they  rely  on  small  growth  rates  and  so  the  unstable  disturbance  of 
concern  must  be  near  to  a  neutral  state.  Thus,  it  was  believed  that  such  theories  are  of  no 
use  in  describing  the  initial  evolution  of  ‘far-from-neutral’  unstable  modes.  However,  several 
recent  studies  have  derived  integro-differential  equations,  using  weakly  nonlinear  theories,  to 
describe  the  nonlinear  evolution  of  (general)  unstable  modes  on  a  variety  of  shear  layers  (see  the 
previous  paragraph).  These  studies  are  based  on  the  assumption  that,  in  actual  physical  flow 
situations,  shear  layer  spreading  or  other  external  changes/effects  would  result  in  the  otherwise 
relatively  unstable  modes  having  their  growth  rates  diminished  in  real  terms,  so  that  a  weakly- 
nonlinear  critical-layer  theory  becomes  appropriate;  in  fact,  we  adopt  the  same  argument  for 
this  present  study.  The  general  proposal/ argument  is  supported  by  the  findings  of  Michalke 
(1964),  Crighton  &  Gaster  (1976)  and  Hultgren  (1992). 

In  this  study,  we  use  the  unsteady  critical-layer  and  weakly-nonlinear  theories  to  describe 
the  nonlinear  development  and  evolution  of  unstable  linear  disturbance  modes  on  stratified 
shear  flows  where  the  Richardson  number  J  takes  values  less  than  one  quarter.  In  particular, 
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we  consider  the  temporal  evolution  of  two-dimensional  (2D)  modes  on  Drazin’s  (1958)  model 
flow;  we  note  however  that  the  analysis  can  easily  be  modified  for  the  spatial  evolution  case 
and/or  a  different  choice  of  model  flow.  In  addition  to  not  complicating  the  analysis,  the  choice 
of  considering  only  2D  disturbances  is  justified  by  noting  (i)  the  model  flow  is  also  only  2D,  and 
(ii),  that  Squire’s  theorem  holds  for  such  flows  where  Boussinesq’s  approximation  is  used  (see 
Koppel,  1964).  Perhaps  the  main  assumption  made,  in  this  first  such  study  for  J  <  1/4,  is  that 
the  Reynolds  number  of  the  flow  is  large  enough  such  that  the  critical  layer  is  not  significantly 
affected  by  viscous  effects  to  the  orders  considered.  Such  an  assumption  appears  justified  based 
on  our  results. 

Although  this  paper  is  solely  concerned  with  the  nonlinear  evolution  of  (non-marginally) 
unstable  stratified  shear  flow  which  has  important  geophysical  applications,  the  actual  initial 
motivation  for  the  study  was  the  authors’  desire  to  develop  a  theory  to  describe  the  nonlinear 
evolution  of  the  inviscid  Gortler  modes  considered  initially  by  Bassom  &  Hall  (1991),  and  lately 
by  Dando  (1992),  Blackaby  &:  Choudhari  (1993).  The  first  two  of  the  above  papers  demonstrate 
that  in  the  presence  of  a  relatively  weak  cross  flow,  longitudinal  vortex  disturbances  of  all 
wavelengths  are  stabilised  such  that  the  inviscid  modes  possess  some  of  the  largest  growth  rates 
whilst  also  being  neutral  at'  certain  other  wavenumbers.  Their  governing  equation  is  similar 
to  the  Taylor-Goldstein  equation  which,  as  mentioned  earlier,  governs  the  linear  stability  of 
stratified  shear  flows.  In  fact  Blackaby  &  Choudhari  (1993)  have  illustrated  the  close  connection 
between  the  two  problems  and  propose  a  definition  of  a  generalised  Richardson  number  for 
such  centrifically-driven  instabilities.  The  ideas  developed  here  in  this  paper  also  have  obvious 
applications  to  inviscid  modes  in  flow  above  a  heated  plate,  similar  to  those  considered  by  Hall 
h  Morris  (1992). 

Whilst  different  from  the  approach  adopted  in  this  study,  there  are  alternate/ complement  ary 
nonlinear  theories  that  have  been  developed  recently  in  which  two  or  more  of  the  flow  distur¬ 
bances  mutually  interact.  Such  theories  generally  require  smaller  disturbance  amplitudes  but 
may  also  need  the  disturbances  to  exist  in  specific  ‘configurations’.  These  other  theories  are 
generally  referred  to  as  wave/wave  and  vortex/wave  interactions.  For  a  discussion  of  wave/wave 
interactions  and  resonant-triads  the  reader  is  directed  to  the  book  by  Craik  (1985)  and  the  pa¬ 
pers  by  Goldstein  &  Lee  (1992),  Wu  (1992).  Strongly  nonlinear  vortex/wave  interactions  were 
first  looked  at  by  Hall  &  Smith  (1991)  and  their  ideas  were  clarified  and  extended  by  Brown, 
Brown,  Smith  &  Timoshin  (1993).  Dando  (1993)  has  looked  at  this  type  of  interaction  in  a 
heated  boundary  layer  where  both  streamwise  vortices  and  inviscid  travelling  waves  are  present. 

The  format  of  the  rest  of  the  paper  is  as  follows.  In  the  next  section  we  present  some 
background  details  of  the  flow  concerning  us  in  this  paper,  namely  the  stratified  shear  flow 
model  due  to  Drazin  (1958),  as  well  as  introducing  some  notation  and  concepts  ready  for  the 
proceeding  sections.  In  §3  the  flow  outside  the  critical  layer  is  considered;  whilst  in  §4  we 
consider  the  flow  inside  the  critical  layer.  In  §5  we  discuss  the  range  of  validity  and  solution 
properties  of  the  evolution  equations.  Finally,  in  §6  we  draw  some  conclusions. 
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2  The  stratified  shear  flow  model 


Following  Churilov  &  Shukhman  (1987,88)  we  choose  to  consider  the  simple  but  realistic 
stratified  shear  flow  model  of  Drazin  (1958),  where  the  unperturbed  velocity  in  the  x  direction, 
u0,  and  the  density,  po,  are  dependent  on  the  vertical  coordinate  y  (in  fact  uo  =  tanh  y  and 
Po  =  Poc  —  V )•  The  governing  equations  of  motion  in  the  Boussinesq  approximation  can  then 
be  written  in  the  form 

—Arp  -  ip)  =  T]KA2ip, 

^  +  { P ,  0}  =  nAp,  (2.1  a,  b) 

where  {a,  6}  =  (da/ dx)(db/ dy)  —  (da/ dy)(db/ dx);  ip  is  the  stream  function;  p  is  the  density;  rj 
the  Prandtl  number;  k  the  thermometric  conductivity  and  J  the  Richardson  number. 

The  shape  of  the  stream  function  of  a  two-dimensional  disturbance  mode  of  the  inviscid 
linear  problem,  denoted  ip[ ^  say,  is  found  to  satisfy  a  Taylor-Goldstein  equation  of  the  form 

Uip[1]  =  0,  (2.2) 


where  the  operators  L i  have  the  form 


_  _  (([,x2  _  2 sinh ysech3y _ J _ \ 

~  dy2  \  (tanh  y-c)  (tanhy-c)2/ 


(2.3) 


Here  k  and  c  are  respectively  the  wavenumber  and  wavespeed  of  the  infinitesimal  distur¬ 
bance,  whilst  the  index  /  corresponds  to  the  harmonic  being  considered  (see  later).  Note  that 
the  L i  are  singular  where  tanh  y  =  c;  at  such  locations,  usually  referred  to  as  critical  levels,  the 
assumptions  employed  in  the  linear  inviscid  approach  are  no  longer  valid  and  special  attention 
must  be  paid.  Thin  regions  around  such  levels,  usually  referred  to  as  critical  layers,  must  be 
introduced  into  the  mathematical  model;  such  layers  axe  generally  the  first  to  ‘feel’  the  effects 
of  increasing  disturbance  amplitude. 

Drazin  (1958)  found  neutral  eigensolutions  of  the  form 

W1*  =  B±tpa(y ),  tpa(y)  =  sinh1-fc*  |y|  sech  y, 

with  c  =  0  and  J  =  fc2(l  —  fc2),  (2.4a  —  d ) 

satisfying  (2.2)  coupled  with  the  boundary  conditions  ip[^  — *■  0  as  y  — *  ±oo.  Here  the  ±  on  B 
are  related  to  vertical  position  with  respect  to  the  critical  level  at  y  =  0  i.e.  y  >  0  and  y  <  0 
respectively.  The  relationship  between  the  amplitudes  B+  and  B-  will  follow  from  matching 
to  the  inner  problem  (cf.  Miles,  1961).  In  fact,  in  §4  we  shall  find  that  f3_  =  i~1~2uB  where 
we  have  written  B  =  B+  and  the  Frobenius  root  v  is  defined  below. 

For  later  convenience,  we  introduce  the  quantity 

v  =  ±^Vl  -4 J; 


4 


(2.5) 


note  that,  unlike  Miles  (1961),  we  choose  to  include  the  “±”  in  our  actual  definition  of  v  — 
this  is  also  done  for  convenience.  It  easily  follows  from  (2.4,5)  that  for  Drazin’s  model  flow 
v  =  1/2  —  k2  and  thus  the  appropriate  root  is  immediately  determined.  However,  for  general 
flows,  an  analytical  solution  is  not  possible  and  the  relevant  choice  of  sign  for  given  k ,  c  and  J 
can  only  be  determined  by  inspecting  the  shape  of  the  eigenfunction  obtained  from  a  numerical 
solution;  however,  there  does  appear  to  be  a  common  pattern  between  location  on  neutral  curve 
and  which  sign  in  (2.5)  is  appropriate  (see  Blackaby  &;  Choudhari,  1993).  Note  that  ‘half’  of 
Drazin’s  neutral  curve  corresponds  to  the  positive  choice  of  sign  and  the  other  ‘half’  to  the 
negative  choice  (see  Figure  2.1);  moreover,  note  that  Drazin’s  solution  is  also  valid  for  J  <  0 
which  corresponds  physically  to  buoyancy  driven  flows.  For  J  <  0,  it  appears  that  the  w— ” 
option  is  always  relevant;  thus  the  eigenfunction  also  is  singular  at  the  critical  level. 

A  study  of  the  fundamental  and  other  low  harmonics  outside  and  inside  the  critical  layer 
is  necessary  to  derive  the  desired  amplitude  equations;  in  the  next  section  we  consider  the 
flow  outside  of  the  critical  layer.  Here  the  details  are  dependent  on  the  flow  model  under 
consideration,  but  the  method  is  quite  general  and  can  be  applied  to  other  flows.  In  §4  we  shall 
see  that  the  critical-layer  analysis  is  almost  entirely  independent  of  the  stratified  shear  flow 
being  considered. 


3  Outside  the  critical  layer 


We  consider  a  wave  of  sufficiently  small  amplitude  in  the  neighbourhood  of  a  general  point 
(J,  k)  on  the  neutral  curve  of  Drazin  (1954)  (see  Figure  2.1).  Note  that,  since  we  are  not 
considering  the  initial  evolution  of  a  mode  on  a  marginally  unstable  flow,  we  are  not  restricted 
to  the  sole  case  J  ~  1/4.  In  fact,  we  begin  by  allowing  J  to  take  any  value  <1/4  (i.e.  we  choose 
to  exclude  the  special  case  J  =  1/4),  and  then  see  which  values  of  v  =  u(J,  k)  require  special 
attention  later.  Note  that  J  /  1/4  for  all  but  one  point  on  each  of  Drazin’s  and  Hazels’  (1972) 
neutral  curves;  moreover,  the  generalised  Richardson  number  is  always  negative  for  neutral 
inviscid  longitudinal  vortices  in  3D  boundary  layers.  In  fact,  the  theory  of  this  and  the  next 
section  can  be  regarded  as  the  extension  of  the  work  of  Churilov  &  Shukhman  (1988)  to  the 
case  v  ^  0,  i.e.  to  a  far  wider  range  of  problems.  At  the  outset,  we  also  choose  to  exclude  the 
special  case  J  =  0,  which  corresponds  to  a  flow  with  no  stratification  (this  case  is  considered 
by  Goldstein  &  Leib,  1988). 

Following  Churilov  &  Shukhman  (1988),  we  introduce  the  small  parameter  c,  characterising 
the  magnitude  of  the  mode.  It  is  also  necessary  to  introduce  the  ‘slow’  evolutionary  time  r  =  fit, 
with  fi  also  small  i.e.  the  amplitudes  appearing  in  (2.4a)  are  considered  to  be  functions  of 
r.  Note  that  although  the  time-scale  r  is  slow  with  respect  to  the  inviscid  timescale  t,  it  is 
still  much  ‘faster’  than  the  timescale  t  ~  k~1  of  the  viscous  spreading  of  the  base  shear  flow, 
provided 

*  <  n  (<  i); 

this  restriction  is  not  severe  due  to  the  large  size  of  the  Reynolds  number  Re  ~  «-1  for  flows  of 
practical  interest.  Later  we  shall  have  to  relate  the  small  parameters  k,  e  and  fi  to  one  another 
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so  that  we  can  balance  the  linear  terms  of  the  amplitude  equation  with  possible  nonlinear  ones. 
We  write 

J  —  Jo  +  pJ\i  (3-1) 

where  Jo  corresponds  to  a  general  point  on  the  neutral  curve  and  J\  <  0  is  order  one.  The 
streamfunction  and  density  are  expanded  as  Fourier  series  in  x , 

OO  OO 

4>  =  In ( cosh t/)  +  ^2  MT>y)e'klx,  P  =  poc  ~  V  +  J2  Pi(T>y)e'klXi  (3.2a, 6) 

/=  — OO  /=  — OO 

with  boundary  conditions 

ipi,  Pi  0  as  y  — +  ±oo.  (3.3) 

In  the  present  study,  we  need  to  consider  the  fundamental,  the  zeroth  harmonic,  the  second 
harmonic  and  the  third  harmonic.  Equations  (2.1a,b)  are  solved  both  outside  and  inside  the 
critical  layer  and  then  these  solutions  are  matched  to  obtain  the  evolution  equations.  Outside 
the  critical  layer  the  fundamental  harmonic  dominates  the  perturbation  and  has  tin  amplitude 
of  order  t.  The  other  harmonics  are  the  result  of  self-interactions,  with  the  zeroth  and  second 
having  amplitudes  of  order  e2  and  the  third  of  order  e3. 


3.1  The  fundamental  harmonic 

In  order  to  derive  the  evolution  equation,  the  following  terms  of  the  expansion 

V>i  =  +  e/iip  j3)  +  e/c^j3)  +  •  •  • ,  (3.4) 

need  to  be  considered.  Here  rj> is  the  neutral  mode  of  the  inviscid  linear  problem;  takes 
into  account  the  r-dependence  of  the  solution;  and  xj>^  is  a  correction  to  tp^  for  dissipative 
(viscous)  effects.  The  corresponding  p\  expansion  includes  analogous  terms.  Thus 

rpP  =  B±(r)iJ>a(v)  where  rpa{y)  =  sinh^"  |y|sechy;  (3.5a, 6) 

we  note  that 

rf1  ~  B±M|sr|i+\  as  y  — *  0.  (3.6) 

The  second  term  in  the  expansion  (3.4)  satisfies 

=  Qu  *P\^  0  as  y  — ►  ±oo,  (3.7a) 


where 


JiB±rpa 


(3.76) 


n  =  _  2 iipa  dB±  (  1  Jo  \ 

1  tanh*  y  fctanhy  dr  \cosh2  y  tanh3  y ) 

The  solution  can  be  considered  to  be  the  sum,  ip[2^  =  °i  a  particular  integral 

of  (3.7),  say,  and  the  complementary  function,  ip[ say-  As  y  — ►  0,  it  follows  from  Taylor 
(1931)  that 

=  B4«SJ2l»|i+v(l  +  0(|!)|-1))  +  B±6S|s,|S-(l  +  0(|s|-')),  (3.8) 
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where  and  b^J.  are  constants  as  yet  undetermined.  Note  that  if  the  Frobenius  roots  |  ±  |i/| 
differ  by  an  integer  then  (3.8)  is  no  longer  appropriate  (logarithms  are  needed).  As  such  cases 
( 1 /  =  |m;  m  integer)  are  isolated,  we  choose  not  to  concern  ourselves  with  them  (and  their 
immediate  neighbourhood)  in  this  paper. 

A  solvability  condition  for  the  above  boundary-value  problem  (3.7)  is  required.  Note  that: 
(i)  the  operator  Li  is  self-adjoint  away  from  the  critical  level  y  =  0,  and  (ii)  the  right-hand  side 
is  singular  at  y  —  0.  Rather  than  closely  follow  the  method  of  Churilov  &  Shukhman  (1987) 
for  deriving  the  modified,  solvability  condition,  for  Jq  ^  1/4  it  is  more  convenient  to  adopt  the 
related  approach  employed  by,  for  instance,  Hickernell  (1984). 

The  solvability  condition  is  derived  by  multiplying  both  sides  of  equation  (3.7a)  by  xpa  and 
integrating  over  all  y,  excluding  the  (sole)  critical  layer  at  y  =  0.  After  integrating  by  parts; 
imposing  the  boundary  conditions  at  y  =  ±oo;  and  the  asymptotic  forms  of  i/>a  and  V’lcV  35 
y  — ►  0,  we  find  that 

jMi  dy  =  -  {'PwF'P'a  -  V’lCF^a]”.  =  (Bib [5  +  Bib[2})  ,  (3.9a) 

where  the  bar  through  the  integral  sign  indicates  that  the  finite  part  of  the  integral  should  be 
taken.  After  substituting  for  Qi  and  using  the  relations  B-  =  i~l~2u B,  B+  =  B  (to  be  derived 
in  the  next  section),  the  solvability  condition  becomes 


4t1+2,,/i 


cos(iri/) 


dB  Jilik  .  .  n 

jr+ir*-0"’)* 


=  -2i/fl(6™  -  r'-b j!>),  (3.96) 


where 


h- 4  8mh,—  ( — ^  dy  and  h  =  4  sinh2l/_1  ydy.  (3.9c, d) 

Jo  coshy  \cosh  y  tanh  y/  Jo 


Four  relations  involving  and  b^2}  are  determined  from  the  inner  problem  considered  in 
the  next  section,  thus  determining  the  possible  evolution  equations  for  the  wave  amplitude 
B{r). 


3.2  The  zeroth,  second  and  third  harmonics 

The  presence  of  these  terms  is  due  to  the  process  of  harmonic  generation,  i.e.  due  to 
nonlinearity  the  fundamental  of  0(c)  generates  the  zeroth  and  second  of  0(e2)  and  so  on.  The 
zeroth  and  second  harmonics  are  expanded  in  the  form 

*  =  cVS”  +  +  ‘Vf  +  •  •  • ,  (3.10) 

*  =  i ty™  +  <VV4J)  +  e2«i(43)  +  •  •  ■ ,  (3.11) 

and  similarly  for  po  and  p2.  It  is  only  necessary  to  consider  the  first  term  of  the  third  harmonic, 

namely 

V»3  =  +  •  •  • ,  (3.12) 
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and  the  similarly  for  p 3.  The  resulting  sets  of  equations  from  expansions  (3.10-12)  do  not  have 
simple  exact  solutions  when  Jo  7^  1/4;  however  it  is,  in  general,  a  relatively  simple  process, 
if  a  little  tedious,  to  deduce  their  asymptotic  forms  as  the  critical  level  is  approached.  The 
significant  terms  of  these  asymptotes  are  quoted  in  the  next  subsection. 

However,  the  first  term  of  the  second  harmonic  requires  extra  attention  as  it  can  lead  to  the 
largest  nonlinearity  in  certain  circumstances  when  v  is  positive  (see  §4).  The  necessary  analysis 
is  simply  the  generalisation  of  that  for  the  case  v  =  0  (J  =  1/4)  presented  by  Churilov  & 
Shukhman  (1987)  in  their  unpublished  Appendix  B.  Since  this  important  work  is  unfortunately 
not  readily  available  for  reference,  we  take  the  liberty  of  presenting,  in  Appendix  A  of  this 
paper,  more  details  of  their  method  than  we  would  otherwise. 

3.3  The  asymptotic  expansions  as  y  — *  0 

In  terms  of  the  new  variable  Y  =  p~xy  and  the  functions 

<j>i  =  fa  -  {1 +  2i/)pi,  (3.13) 

where  pi  =  —2 dpi/dy,  the  asymptotes  for  the  fundamental,  second,  zeroth  and  third  harmonics 
as  y  — ►  0  are 


fa  =  e^+' 


2k 


’+*'  I a^B±  - 


3t  Y 

J\B± 


+  t^-''b(?lB±\Y\ 


2 1/ 


In  1^1]  m 


4+*- 


+e'C"  *  48* 


{*7(5  —  2i/)  +  (1  —  2i/)}B±~^ - \- 


Y 3 


fa  =  ejiH*'.  0  + 


-4 *g*6j± 

(1-2 v)  1  1 

B  « 


+  e^+v 


-JiB± 

v{\+2v) 


+€K(l 


-|+» 


*(*?  -  1)(1  -  4v2)B±  |yr|i+*1 


8k 


Y3 


(3.14a,  6) 


fa  =  C  V1+3* 


(1  +  2 v)B\  \Y\1+iv' 
2(3  -  2i /)  Y 2 


+  •  •  ■ ■  +  +  •  •  • 


+t2fi^Blg2±\Y\^  +  t2^-vBlh2±\Y\^  +  •  •  • , 


8 


<f> 2  =  e*fj~1+3l/. 0  -  ein 


2..2v 


JiBl  \Y\1+2l/ 
2z/(3  -  2v)  Y2 


+  •  ■  ■  -  „h- 


4i/c+ 


-e2n^Blh 


Y 5 


4  v 


(1-2*) 
\Y\l-*  +  - 


\YV 


(3.15a, b) 


dtpp 

dr 


cV1^" 


(1  +  2v)  d(l?±|2  |K|1+21' 


dr  Y2 


Wo  -  ,2  -i+2,  n  _  .2s±2»Ji  d\B±\2  |*T+2* 

0r  ^  ^  2*  dr  Y2 


+e2KH~*+2l/ 


|(1-4k2)(1-2k)(,-1)|B*|3£1j 


l+2i/ 


+ 


(3.16a, 6) 


_  ,  tw>(H-2o)(l-2y)  Bj 

™  H  2(3  + 2k)3  |K|!-3"  +  ’ 


*  =  <V*+3‘'-  <>+•••  , 

(3. 17a, 6) 

Note  that  above,  for  the  sake  of  brevity,  we  have  only  included  the  more  important  and 
illustrative  terms  with  regard  to  our  aim  of  deriving  amplitude-equations  for  B.  It  should  be 
noted  that  the  above  asymptotes,  in  general,  only  contain  the  leading  term  (for  Y  —*  oo)  at  each 
order.  Further,  it  should  also  be  noted  that  these  (incomplete)  asymptotes  may,  depending  on 
the  relative  sizes  of  e,  k  and  /x,  be  disordered  as  written;  it  is  sensible  to  postpone  a  discussion 
concerning  the  relative  sizes  of  these  quantities  until  §5,  when  the  sizes  of  the  competing  nonlin¬ 
earities  are  known.  Here  <f±,c±  and  q±  are  constants,  as  introduced  by  Churilov  &  Shukhman 
(1988);  the  first  two  are  in  fact  expansions,  each  term  of  an  order  yet  to  be  determined.  The 
constants  g?±  and  fi2±  are  defined  in  Appendix  A  of  this  paper. 

It  is  worthwhile  to  consider  these  asymptotes  a  little  further.  We  see  that  the  leading  terms 
in  the  ^-expansions  are  all  zero;  in  fact,  such  behavior  is  in  full  agreement  with  the  theorem 
of  Miles  (1961).  In  the  next  section,  we  shall  see  that  this  fact  results  in  the  ‘biggest’  cubic- 
nonlinearity  not  contributing  to  the  amplitude  equation;  the  same  result  is  also  found  when  J  = 
1/4  {v  —  0)  and  is  responsible  for  the  quintic  nonlinearities  occurring  in  the  evolution  equations 
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derived  by  Brown  &:  Stewartson  (1978),  Churilov  &  Shukhman  (1988)  and  in  this  paper.  As 
pointed  out  by  these  previous  authors,  for  a  ‘contributing’ jump,  it  is  in  general  necessary  for  the 
(^-asymptotes  to  be  non-zero.  They  note  that  this  condition  is  satisfied  if  the  complementary- 
function  terms  in  the  <f> 2  asymptote  are  non-zero  (thus  resulting  in  a  ‘contributing’  quintic 
nonlinearity).  Churilov  &;  Shukhman  (1988)  also  note  that  if  viscosity  effects  are  included  in 
in  the  original  formulation  of  the  problem,  then  the  above  condition  is  satisfied  as  long  as  the 
Prandtl  number  77  ^  1  (thus  resulting  in  a  ‘contributing’  cubic  nonlinearity).  In  addition  to 
these  two  so-called  symmetry  breaking  mechanisms,  it  is  clearly  obvious  from  the  third  term 
of  (3.14b)  that  we  have  an  additional  mechanism  for  non-zero  1/,  namely  non-zero  J\.  This 
novel  mechanism  is  made  possible  by  the  two  Frobenius  roots,  1/2  ±  \v\,  not  differing  by  an 
integer.  As  far  as  the  authors  are  aware,  such  a  nonlinear,  non-equilibrium  critical-layer  has 
not  been  the  subject  of  any  previous  studies.  This  ‘new’  mechanism  is  essentially  responsible 
for  the  ‘non-viscous’  cubic  jump  term  in  our  later  amplitude  equation  when  v  is  negative  (later 
referred  to  as  the  ‘Jj-cubic’  term). 

The  fourth  symmetry-breaking  mechanism  possible  here  corresponds  to  that  cleverly  identi¬ 
fied  by  Churilov  &  Shukhman  (1987)  in  the  unpublished  Appendix  B  of  their  paper.  Ironically, 
there  is  no  mention  of  it  in  the  sequel  paper  (Churilov  &  Shukhman,  1988);  however,  it  is  of 
lesser  importance  for  disturbances  of  marginally  unstable  flows  (see  later  for  an  extended  dis¬ 
cussion  Dn  where  each  of  the  four  possible  nonlinearities  in  the  evolution  equation  for  B  are  the 
most  significant).  The  important  thing  to  note  at  this  stage  is  that  the  term  involving  /i2±  in 
the  <j) 2  asymptote  is  not  at  one  of  the  orders  to  be  considered,  in  §4,  for  the  preceding  term  in¬ 
volving  c2±.  In  the  next  section,  we  shall  see  that  this  symmetry-breaking  mechanism  provides 
the  ‘non-viscous’  cubic  jump  when  v  is  positive  (later  referred  to  as  the  outer-complementary- 
function  [OCF]  cubic  term). 


4  The  Critical  Layer 

As  usual  in  such  nonlinear  studies,  the  main  purpose  of  this  section  is  to  calculate  the 
second  relations  between  b\2+  and  b^l  (the  first  being  given  by  the  solvability  condition,  (3.9b)) 
and  thereby  obtain  the  desired  nonlinear  evolution  equation(s)  for  the  disturbance  amplitude 
B(t).  Following  Brown  &  Stewartson  (1978)  and  Churilov  &  Shukhman  (1988),  we  define  new 
functions  It,  P  and  $  where 

V>  =  5 fi2Y2  +  V,  p  =  p0-fiY-2fx-1PY ,  $  =  -  (1  +  2v)P.  (4.1a  -  c) 

Inserting  these  into  the  governing  equations  (2.1),  (2.2)  gives 

=  1(1  -  2*0*.  -  *}*  "  ~  *«)  +  W3*yyy,  (4.2 a) 

%.„)*  =  -/T2{*y,¥}*  -  ~  *.)  +  *H~3*yyy  +  (v  ~  1  )kh~3*yyy,  (4.2 b) 
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(4.3) 


where  {a,  b}*  =  axby  —  ay&x >  the  operator 


and  it  is  assumed  that  Kfi~3  <C  1  i.e.  viscous  effects  are  not  large  enough  to  affect  the  operator 
Nx  at  leading  order. 

In  the  rest  of  this  section  we  proceed  to  solve  equations  (4.2)  for  the  relevant  lower  order 
terms  of  the  lower  harmonics.  In  fact,  we  need  to  consider  the  cases  v  >  0  and  v  <  0  separately; 
this  is  because  the  largest  ‘non-viscous-in-origin’  cubic  term  in  the  corresponding  evolution 
equations  is  different  in  these  cases.  When  v  <  0,  this  term  is  proportional  to  J\\  but  when 
v  >  0,  this  term  is  related  to  the  complementary  function  of  the  second  harmonic  of  the  inviscid, 
outer  problem.  Note  that  when  v  >  0,  the  second  and  third  terms  of  the  ipi  asymptote  (3.14a) 
are  ordered  and  thus  the  cubic  nonlinearity  proportional  to  J\  would  be  insignificant  compared 
to  that  proportional  to  b^±-  Note  that  the  other  two  significant  nonlinearities  (a  cubic,  due  to 
viscosity,  and  a  quintic,  due  to  complementary  function  terms  of  the  critical  layer  solution  for 
the  second  harmonic)  are  possible  in  either  case. 


(i)  The  case  v  <  0 

The  solution  is  again  constructed  in  the  form  of  a  Fourier  series  in  z;  we  expand  the 
fundamental,  zeroth,  second  and  third  harmonics,  respectively,  as  follows 

=  e/i^+,,^i1)  +  •  •  •  +  J2)  +  •  •  •  +  e/<a'H'tf,(13a)  +  •  •  •  +  +  •  •  • 

+eV"!+3,/^S4)  +  •  •  •  +  CK/i-3+*'^(1s)  +  •  •  •  +  e3Kjr^+3,/$S6)  +  •  •  •  +  t*n~ ^+5,/$S7)  +  •  •  • , 

=  e2fi~1+2u¥o)  +  •  •  •  +  £V*'tf£,)  +  *  •  •  +  e2KH~4+2,,'I/{o)  +  •  •  •  +  e4/i'4+4*'^7)  +  •  •  • , 


=  e2  +  •  •  •  +  23)  +  ■  •  •  +  e2K/r4+2"'If23)  +  •  •  •  e V"4+4,/$i7)  +  •  •  • , 


*3  =  eV*+3*'^1)  +  •  •  • , 

(4.4a  —  d) 

and  similarly  for  the  $,’s.  These  expansions  are  not  necessarily  completely  ordered  (depending 
on  the  sizes  of  £,  k  and  /i);  moreover,  only  those  terms  crucial  to  deriving  the  evolution  equations 
have  been  included.  The  scalings  follow  directly  from  the  outer  asymptotes  (3.14-17)  and/or 
by  considering  the  process  of  harmonic  generation. 


4.1  0(e/i5+,/)  of  the  fundamental 


Upon  writing 


we  find  that  at  this  order  equations  (4.2a,b)  give 


(4.5) 


=  4(1- S-)**!1’,  =  0  • 


(4.6a,  6 ) 


A  solution  of  these  equations  which  matches  to  the  outer  solution  (3.5)  is 

(i  +  2i/)r(i  +  i/)i*-3‘' 


^  =  W(r,Y)  = 


/  t  »  vB{t  —  t)e 
Jc 


—ikYt 


dt, 


(4.7a) 

(4.76) 


Airkk+1'  Jc 

*1"  =  0, 

where  the  contour  C  is  shown  in  Figure  4.1a.  The  function  W(r,Y)  has  a  single  asymptotic 
representation  in  the  lower  half-plane  (— ir  <  aigY  <  0) 

W(T,Y)  =  B(T)Yt+‘'  +  0(V->+,/),  as  \Y\  oo.  (4.8) 


Later  we  derive  evolution  equations  for  the  amplitude  B(t)  but  for  the  moment  we  can  regard 
it  as  an  arbitrary  function  that  satisfies  the  requirement  B(r)  — »  0  as  r  — »  — oo.  Matching  with 
the  ‘outer’  asymptote  (3.14a)  yields 

B-  =  (-B+  =  B);  (4.9) 


this  result  has  already  been  referred  to  and  used  in  the  previous  sections. 


4.2  0(e2 fi~1+2v)  of  the  zeroth  harmonic 

At  this  order  equations  (4.2a, b)  yield 

~^{oV  =  ik{wW  -  WW')\  =  0,  (4.10a, b ) 

dr  or 

where  the  overbar  again  denotes  the  complex  conjugate  and  prime  denotes  differentiation  with 
respect  to  Y.  A  solution  of  the  above  which  matches  to  the  outer  solutions  is 

*i>"  =  (1  |2^|H,T.  *?>  =  o.  (4.11a, 6) 
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4.3  O (e2/i  1+2")  of  the  second  harmonic 


At  this  order  equations  (4.2a, b)  yield 

Nttf+r) =  ik(  1  -  2i/)$<1}  +  ik(WW"  -  W'2)  and  =  0.  (4.12a,  6) 

We  note  that,  as  the  right-hand  sides  of  (4.12a,b)  do  not  involve  W,  and  have 
unique  asymptotic  representations  as  |  V|  — ►  oo  (in  the  lower  half  plane  of  complex  Y )  of  the 
form 


~  ^±^1  B*Y-w  +  (m^Y^  +  n^Y^)  +  •  •  • , 


Ai/m^ 

(T^) 


y$-L 


+ 


(4.13a,  6) 

Matching  with  the  asymptotes  (3.15a, b)  requires  that  c±  =  eV_t+3,/ci±,  so  that  Ci+  =  m^\ 
Ci-  =  t*1+2l/mj1),  n2l,c1+q+  =  and  fi^c^q.  =  r1_2*'nj1).  These  lead  to  the  relation 

Ci_(?-  —  i~4t/q+)  =  0.  However,  as  the  adopted  model  flow  is  symmetric,  q+  =  and  so,  for 
v  ^  0,  Ci_  =  0  =  =  n^\  Thus 


(4.14a) 


and  the  solution  to  (4.12a)  is 


B(T-t-t2)(ht2)-$-‘'(t1-t2)2(ti+t2)$+,'(2t+tl+h)-i-,'e-ikY{2t+t'+h\  (4.146) 


4.4  0(eV  f+3,/)  of  the  fundamental 

Prom  (4.2a, b), 

=  *^(1  “  2«/)$(,2)  +  ik(9iirW  -  iyW’  +  2*PW"  -  ¥?'W  -  WW), 

^M^)*i2)  =  0.  (4.15a,  6) 

Due  to  the  earlier  $  terms  being  zero,  the  right-hand  side  of  (4.15b)  is  also  zero.  Thus 
eJ2V,  where  ej2)  is  a  constant  and  V  satisfies  IV, ,(£_„)  V  =  0  with  the  asymptotic  representation 
V  ~  Y*~v  as  |K|  ->  oo.  Any  non-zero  choice  of  e[2)  in  (4.15a)  results  in  a  Yi~v  term  in  the 
asymptotic  representation  of  $(,2)  and  this  contradicts  the  theorem  of  Miles  (1961),  from  which 
we  know  that  the  expansion  of  ij)a  as  y  — ►  0  cannot  involve  a  term  proportional  to  y^~v .  Thus 
it  follows  that  $,2^  =  0. 
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We  also  note  that  the  right-hand  side  of  equation  (4.15a)  contains  W  and  so  'I')  '  does 
not  have  a  unique  asymptotic  representation  as  jy|  — >  oo.  However,  because  the  asymptotic 
representation  of  does  not  contain  a  term  proportional  to  it  is  impossible  to  obtain 

non-zero  5^  at  this  order.  Thus  the  strongest  nonlinearity  does  not  contribute  to  the  evolution 
equation. 

The  solution  for  '5^  is  needed  for  later  calculations  regarding  the  quintic  nonlinearity;  the 
solution  of  (4.15a)  is 


*F  =  r3q+")  rit  /  itl  /  ih  i  ih 

32  7T3  Jo  Jc  Jc  Jc 

B{t  —  t  —  t\)B{r  -  t  -  U)B(t  -  t  -  l3)(<lh<3)"’"'(<3  -  *1  -  t2)±+u(t3  -t-t !  -  t2)->-l'x 
f<3  ^l(<l  ~  t3)2  +  ^2^2  ~  ^3)2  +  —  t\  —  ^)j  —  *i+2l/^2(tl  +  ^2)  2  V  X 

(*l-*2)2-"(*l+<2+*3)(2t3-*l—  <2)T’^---I/,  ^f))/f(ti  -t2) 

(4.16) 

Here  H(x)  is  the  so-called  Heaviside  function: 

ffM=U::x<°0, 

and  F(a,  6;  c;  z)  is  the  hypergeometric  function  (see  Erdelyi,  1953;  Abramowitz  &  Stegun,  1964). 


4.5  0(e/i§+1/)  of  the  fundamental 


A  major  difference  of  our  problem  for  general  v,  compared  to  the  v  =  0  case  considered  by 
Churilov  &  Shukhman  (1988),  is  that  there  is  no  linear  contribution  to  the  evolution  equation 
from  the  critical  layer.  This  can  be  seen  from  the  analysis  of  this  and  the  next  subsection. 

At  order  t^+v  terms  proportional  to  Ji  first  enter;  equations  (4.2a,b)  yield 


which  have  solutions 


*uj+„)*(.3-)  =  j(i  -  2»)i»l3,)  -  W. 

=  —  (1  +  2u)W' 


(4.17a, 6) 


(4.18a, 6) 


Note  that  for  definiteness,  we  choose  to  include  the  complementary  function  alongside  the 
nonlinear  term  which  is  going  to  be  balanced  with  at  this  order. 
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4.6  0(efi2~,/)  of  the  fundamental 

Here  we  find  that  (4.2a, b)  yield 

=  j(l  -  Nuh^b'  =  0.  (4.190,6) 

These  have  solution 

6)  =  0,  $S36)  =  0,  (4.20a,  b) 

where  once  again  we  shall  include  the  associated  complementary  function  alongside  the  non¬ 
linear  term  which  is  going  to  be  balanced  with  at  this  order.  The  solutions  (4.20)  imply  that 
the  whole  of  the  linear  contribution  to  the  evolution  equation  comes  from  outside  of  the  crit¬ 
ical  layer;  note  that  the  logarithm  occurs  in  (4.18a)  and  not  in  (4.20a)  i.e.  any  jump  induced 
by  the  presence  of  the  logarithm  does  not  occur  at  the  desired  order,  O (e/z"*'),  necessary  to 
affect  the  relationship  between  6^  and  6^.  The  authors  have  verified  that  this  is  in  fact  the 
case  by  computing  ‘near-neutral’  linear  growth  rates  from  a  numerical  solution  of  the  Taylor 
Goldstein  equation  (2.2)  and  comparing  them  with  those  predicted  analytically  by  equating  the 
right-hand  side  of  expression  (3.9b)  to  zero. 


4.7  0(  e2H2v)  of  the  zeroth  harmonic 

At  this  order  we  only  need  to  calculate  (4.2b)  gives 

d  ,(2)<  ikJi  nirTi7>  jj 


with  solution 


(4.21) 


(4.22) 


4.8  0(  e2/r")  of  the  second  harmonic 

Equations  (4.2a, b)  yield 

=  **(1  -  2 i/)$(22)  -I-  ik(W"^a)  -  2VK'$S3o)'  +  W^S3o)")  - 


MkJ\ 


3-<*+,/>  3  v  '  3  ’  v  1  11  '  (1  +  2 *0  3  ’ 
N^i-^  -  -^nk)(ww"  -  w'w'>  ~  (TP. sj*?’-  (4-23“'  6> 

For  further  calculations  only  is  needed;  using  the  identity  N2^i+v)  =  N2^ !_„)  —  4ifci/  and 
comparing  the  right-hand  sides  of  (4.12a)  and  (4.23b)  it  is  easy  to  deduce  that 

■  (4.24) 


1/(1  +  2v) 


(4.24) 
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4.9  O  (e3/t  §+3,/)  of  the  fundamental 


At  this  order  the  cubic-jump  proportional  to  J\  emerges.  The  governing  equations  yield  in 
concise  form 

JW.*?’  =  5(1  -  2")i^S41  -  +  «‘.41. 


•d+1/)  1  ~  2^  (1  +  2v) 

N  ,  _  2ikJl  +  o(4) 


(4.25a, fe) 


(1  +  2i/) 

where,  in  particular, 

44)  =  *jb(?53o)'«51)'  +  2$<13a)"®?)  -  +  ^"w  -  $42)"W  -  2$^2)'W') 

=  lkJl  (-2W"vP  +  ^'^iy  -  +  ^'iF  )  .  (4.26) 

“f*  lv) 

To  derive  an  evolution  equation  the  asymptotic  representation  of  4^  is  needed.  This  is 

tf<4)  =  +  D{£)Y±-v  +  0{Y-^U)  as  Y  — ►  ±00  ,  (4.27) 

where,  in  particular, 

D(+]  -  D(y 


ii-vk^-v  r°° 


j— -  r  dY  r  dt  <i-"44)(7-  -  Y)e-,i:y<. 

k  —  l/)  J-o o  Vo 


2^r(j  -  *) 

i(2) 


Note  that  the  term  proportional  to  '&i2)  in  (4.25b)  only  contributes  to  the  jump  C+ *  —  ClY. 
After  substituting  for  and  some  calculation,  we  find  that 

n(4)  r»(4)  i1~4,'k3~4,'Ji(l+2u)T3(~  +  u)  f  f  f  ^  u  ^ 

C+  “C-  = - swfil  -  „)  - Jc*>Jc*>JcMUUt3)  x 


(4.28) 


*«) 


( 


32*2i/2r(±  -  1/) 

fi(r  +  <1  —  t$)B(r  +  <2  —  tz)B(r  4-  t\  +  <2  —  ^3)i2^3{^3  ~  ^2)^3  +  ti  —  ^2)* 
(^3  —  ti  —  <2)*  v  H(t$  —  ti  —  <2) 

+- — - — -  f  dtB{r  +  ti  +  t  —  tz)B{r  +  + 1  —  tz)B(r  +  t\  +  <2  +  2<  —  2t^)x 

4  ./o 

(ti  —  ^2 )2(^  1  ■+■  t2)  =  +l,(2t  T  ti  +  ^2)  2  ^(2^3  —  2 1  —  t\  —  ti)x 
(2t  +  ^1+^2  +  ^3)^3  —  2t  —  ti  —  <2)*  v H{tz  —  2t  —  t\  —  <2)^ 


(4.29) 


Matching  the  ‘inner’  asymptote  (4.27)  with  the  ‘outer’  asymptote  (3.14a)  leads  to  the  rela¬ 
tions 

D{+]  =  b$B+,  DL4)  =  il~2vb(?]B- , 

which  combine  to  give 

D(i]  -  D(_4)  =  £(&$  -  mi2])* ,  (4.30) 

where  the  subscript  Jj  denotes  the  cubic  part  of  the  total  jump  which  is  proportional  to  J\ . 
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4.10  0(c3^  f+3*/)  of  the  third  harmonic 

At  this  order  the  governing  equations  yield 

AT3, (!+„)’ *iX)  =  ^(1  +  2*/)$<1)  +  ik(vP"w  -  m£YW'  +  2*W\V"), 

=  0*  (4.31a,  6) 

The  right-hand  sides  do  not  contain  W  and  so  have  a  unique  asymptotic  representation.  Match¬ 
ing  with  the  solution  outside  the  critical  layer  gives  =  0,  whilst  solving  (4.30a)  gives 


(1)  lu>fc*  3l/(l  +  2v)2r2(\  +  v)  r°°  r  r 

*3  i6x2ra-i/)  J0  dtJ0  dtl Jcdh JcdhB(  1  *x)' 


B(r-t-  t2)B(r  -  t  -  tz)t\^t2  ’  %  *  "(<!  +  t2)-i-‘/(f1  -  x 
(*i  +  <a  +  t3)^+l/(3*  +  <1  +  t2  4-  *3)  ^_l,(ti  +  <2  —  ^3)(^i  + 12  —  2 13)  x 


*/  3 
2’  4 


-ikYQt+U+ti+t^ff^  _ 


(4.31c) 


4.11  O  (e4fi  4+4l/)  of  the  zeroth  harmonic 

At  this  order  it  is  only  necessary  to  determine  Equation  (4.2b)  gives 

■fc*oY  =  0,  (4-32) 

/»W 

which  has  the  solution  $0  ’  =  f(Y).  However,  choosing  a  non-zero  /(F)  is  merely  equivalent 
to  taking  a  different  density  profile  for  the  unperturbed  flow  and  thus,  without  any  loss  of 
generality,  we  choose  to  take  f(Y)  =  0  and  =  0. 


4.12  O (c4/i  4+4")  of  the  second  harmonic 

At  this  order  the  equations  (4.2a, b)  can  be  written  in  the  form 

"  **(1  -  2j/)$(37)  =  ikRs  =  ik(Rn  +  Rm  +  Rz  1), 

=  0,  (4.33a,  6) 
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where 

Rn  =  (WV  S2)"  -  2  W'¥iY  +  V^>S2)),  R02  =  2(*J1)"*?)  -  *J1)'*?)') 

and  R31  =  (SW"**1*  -  2F'^|1)'  -  (4.34a  -  c) 

Following  Churilov  k  Shukhman  (1988),  is  written  as  the  sum  =  ^2 In  +  ^2 l  °f 
a  particular  solution  (iV2a+„)'l!'^  =  ikRs)  and  a  general  solution,  °f  the  associated 

homogeneous  equation.  Denning 


we  find  that 


$<7)  =  /*-2,/coo(l/2  -  1/, 2 k)  [  t~2^C(T  -  t)e~2iktYdt , 

Jc 

=  coo(l/2  +  ^,2fc)  J  -  t)e~2iktYdt 

2-)coo(l/2-v,2fc)  /  r  §+"<7(t  -  t)e~2iktYdt,  (4.35  a, 6) 

where  C(r)  and  D(r)  are  arbitrary  functions  which  tend  to  zero  as  r  -»  —00,  and  the  /i-2" 
factors  have  been  added  for  later  convenience.  As  Y  —*  ±00, 


*£}  =  A/±y»+*'  +  o(y-5+"), 


where 


afJm Z,)£.*rr***  - t’Y^'-vi,r  * »• 


(4.36) 


J  '  (i  +  2i')r(|  +  v)  Jo  •  '  '  '  ' 

Thus  the  asymptotic  representations  of  an^  $2  have  the  forms 

^7)  =  [D( t)  +  M±(t)]Y -  ^1^1c{t)YL2-v  +  ■  •  • , 

$<7)  =  fi-2‘'C(T)Yl-‘' + 

(4.37  a,6) 

To  match  with  the  solution  outside  the  critical  layer,  in  (3.15)  it  is  necessary  to  set  c±  = 
e4^~l+3,/C2±,  leading  to  the  relations 

c2+  =  i'~2vC2-  =  C(t),  c2+q+  =  D(t)  +  M+(r)  and  c2-q~  =  il+2v[D(r )  +  Af_(r)]. 

Upon  setting  q+  =  q  =  q-,  we  find  that 


C(r)  =  -^-csc(y)(M+-M_), 


(4.38) 


and,  as  M+  —  M_  ^  0,  this  fully  determines  $27*  ^  0.  Thus,  the  symmetry  $  =  0  has  been 
broken  and  a  ‘jump’  will  occur  at  the  next  order  in  which  appears  in  the  right-hand  side 
of  the  equations.  The  explicit  form  of  M+  —  A/_  is  given  in  Appendix  B  of  this  paper. 
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4.13  0(e5/i"V+5,/)  of  the  fundamental 


(4.39a,  b) 


At  this  order,  the  governing  equations  are 

«i7) = -  2,)*r> + m7), 

=  -ikfi21'  ($57)"^(i1)  +  2^7)'tf(,1)')  = 

The  asymptotic  representation  of  has  the  form 

«i7)  =  +  d(1)y^-v  +  o(y-»+,/), 

where,  in  particular, 

^  -  d-’ = l  iv  r  duh^iT  -  ‘’v)' 

From  (4.7a)  and  (4.35a), 

Jtf»  =  )“»(*’)  /  dh  I  -  ll)e-^(2,-,)C( r  -  U)-E(r  -  f,), 

2iir  Jc  Jc 

(4.41) 

and  simple  manipulations  give 


-ikYt 


(4.40) 


2if-3"jt!-l'(l-2i/)7rf 
I'MVl 


D?  -  D™  =c,j r  r  -  t)B(r  -  2 t)it,  c,  =  , 


r(|  -  x)r(i  +  f)r(i  -  *) 


(4.42a,  6) 


The  quintic  contribution  to  is  found  by  matching  the  above  asymptote  with  the  outside 
expansion  of  the  fundamental  near  the  critical  layer.  This  process  yields 


D{1]  -  D(_7)  =  B(b{?l  -  i~Avb^l)q, 


(4.43) 


where  the  subscript  q  denotes  the  ‘quintic’  part  of  the  toted  jump.  From  (4. 38), (4.42)  and  (4.43) 
we  obtain 


B(l&  -  =  -!^1cSc(y)j(”a+'[M+(r-()-M.(r-()]B(r-21)<«.  (4.44) 


2? 
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4.14  The  cubic  nonlinearity  due  to  viscosity 

The  remaining  terms  in  the  expansions  (4.4),  proportional  to  /c,  are  driven  by  dissipative 
(viscous)  effects.  They  are  not  of  major  concern  in  this  paper  but  are  included  for  completeness; 
the  equations  which  they  satisfy,  their  solutions  and  further  analysis  can  be  found  in  Dando 
(1993)  (see  also  Churilov  &  Shukhman,  1988).  Providing  the  Prandtl  number  is  not  unity,  a 
non-zero  jump,  —  Z>L6^]o  +  —  D^] 2  say,  arises  at  0(e3/c^~ t+3*')  where 

B{t  —  t\  +  t2)(t\t2t3)  2  11  {t  1  —  <2  —  la)2  U{tl  —  ^2)2  |^2  JQ  dt^B{T  —  2t\  t2  -{■  t3  —  t\)x 
B(t  —  fi  +  <3  —  <4)  ^(1  +  2i/)(t*  4- 1\)  —  2(3  —  2v)t\t2(t2  —  ti)2^ 

—B{t  —  2t\  +  ti  +  t3)B(r  —  t\  +  t3)t\t2  ^(<1  —  t2)3  —  2l3(li  —  t2 )2  +  2l3^  j  H(t\  —  t2  —  I3), 


rn«_n«i  +  + 

1  +  '  j2  256i/7T2(3  -  2v)ra  -  u) 


f  dt\  f  dt2  f  dt3  /  dt 
Jc  Jc  Jc  Jo 


B(t  —  2t\  +  t2  +  t3  +  2 t)B(r  —  t\  +  t3  +  t)B(r  —  t\  +  t2  +  I)(lit2l3)  2  v  x 
(li  —  t2  — ■  <3  —  21)5  "(<3  —  t2)3(t2  +  t3)^^l'(2t  +  I2  +  I3)  *  ^(2^i  —  t2  —  t3  —  21)  x 

^2/j  +  (21  + 12  +  ^3)1+2l/(^2  —  ^)2(^2  +  ta,)1~2t'  +  (21  + 12  +  ^3)*^ H (t\  —  t2  — 13  —  21).  (4.45a,  6) 

Matching  with  the  outer  inviscid  solution  requires  that 

Df  -  T>L6)  =  B{b(?l  -  i~Avb(3l)vy  (4.46) 

where  the  subscript  v  denotes  the  cubic  part  of  the  total  jump  which  is  due  to  dissipative  effects. 

It  is  worthwhile  to  consider  the  effects  of  viscosity  a  little  further  at  this  point.  At  the  start 
of  this  section,  it  was  pointed  out  that  we  assume  «  <C  p3,  so  that  the  effects  of  viscosity  do 
not  enter  the  crucial  critical-layer  operator  Nx  at  leading  order.  Instead,  the  effects  of  viscosity 
occur  as  inhomogeneities  at  certain  lower  orders  of  the  hierarchy  of  critical-layer  equations. 
This  approach  is  entirely  rational  as  long  as  /c//i3  <  1.  However,  to  derive  amplitude  equations 
valid  over  a  larger  range  of  k  values,  one  must  follow  the  approach,  introduced  by  Haberman 
(1972),  of  introducing  a  new  parameter,  \n  say,  where 


Atf  =  K/fi3, 


(4.47) 
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is  taken  to  be  order  one  during  the  critical-layer  analysis.  The  operator  Nx  would  need  modi¬ 
fying: 

Nx  *  Nx  ~  tI^h'qy3  (4.48) 

and  hence  all  the  critical-layer  solutions  would  be  modified.  We  have  chosen  not  to  adopt  this 
approach  here  as  it  would  complicate  the  presentation  of  our  theory  and  analysis;  whilst  it  is 
not  really  necessary  so  long  as  A h  is  reasonably  small  (see  Goldstein  &  Leib,  1989). 


(ii)  The  case  u  >  0 


Clearly  the  viscous-cubic  and  the  quintic  nonlinearities  are  still  possible  when  v  is  positive. 
Moreover,  the  analysis  and  the  evolution  equations  are  exactly  the  same  as  for  the  negative  v 
case.  However,  when  v  is  positive  >  e/is+*'  and  the  ordering  of  terms  in  expansion  (4.4a) 

strictly  needs  to  be  changed  so  that  the  term  is  lower  order  than  the  term  It  is  now 

no  longer  rational  to  balance  the  cubic  term  formed  from  the  4^3q)  term,  with  the  'I'j3fc)  term. 
Thus,  the  Jj -cubic  nonlinearity  considered  for  negative  v  is  no  longer  a  possibility. 

There  are  two  other  possibilities  that  we  can  consider  for  positive  v.  Firstly  there  is  the 
cubic  formed  by  the  'fj3  ^  term,  balancing  with  the  \I^36*  term.  However,  this  among  other 
things,  necessarily  leads  to  the  fully  nonlinear  critical  layer  problem.  Alternately,  there  is 
another  cubic  nonlinearity  (  referred  to  here  as  the  outer-complementary-function  [OCF]  cubic 
term)  similar  to  that  considered  by  Churilov  &  Shukhman  (1987)  (for  unity  Prandtl  number) 
in  the  unpublished  Appendix  B  of  their  paper.  This  cubic  nonlinearity  was  not  considered  for 
negative  v  as  the  Ji -cubic  nonlinearity  is  always  larger.  It  arises  from  considering  the  part  of 
the  complementary  function  term  in  the  asymptote  of  the  second  harmonic  outside  the  critical 
layer  (equation  (3.15a))  to  be  at  an  order  fixed  by  the  outside,  inviscid  problem;  rather  than 
just  considering  it  at  orders  fixed  by  the  process  of  harmonic  generation  inside  the  critical  layer. 

It  is  necessary  to  consider  two  additional  terms  in  the  (slightly  re-ordered)  expansions 
(4.4),  namely  terms  the  expansions  of  the  second  harmonic  and  terms 

4>j8^)  in  the  expansion  of  the  fundamental. 


4.15  0{e2n*~v)  of  the  second  harmonic 

At  this  order  equations  (4.2a, b)  yield 

Nu =  (1  -  2 v)ik$W  +  ik(WW"  -  W'2),  and  Nux_v)^  =  0.  (4.49a  -  b) 

The  solution  for  is  written 

^8)  =  Coo(l/2  -  i/, 2 k)  Jct~ 2+^ n2(T  -  t)e~2iktYdt,  (4.50) 
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where  the  function  n2(r)  is  determined  below.  For  large  |K|, 

~  na(r)Ki-^,  tf<8)  ~  m2(r)yi+*'  -  i1"- Mn2(T)Y$-',  (4.51) 

\v 

where  m2(r)  is  a  suitable  arbitrary  function.  Matching  with  the  expansions  outside  of  the 
critical  layer  yields 


<72-  =  *1+2i/02+,  h2-  =  i1+6l,h2+  and  B2h2+  =  -  ---- ^n2. 

4i/ 

From  the  relations  (4.52a-c)  and  (A.5d),  it  follows  that 

,  .  _  4iH1'i/(a/?a  -  /3cxg)  sin[ff(l  +  2t/)/4]  2  . 

^2\TJ  / 1  a  \  •  /  \  (T), 

(1  -  2i/)aa  sin(7ri/) 

and  thus  as  is  non-zero,  it  will  lead  to  a  ‘contributing’  nonlinear  jump. 


(4.52a  -  c) 


(4.53) 


4.16  0(ezfi~l)  of  the  fundamental 

The  governing  equations  at  this  order  are 

Ni,^+u)^i)  =  ^(1  ~  2t/)t*$i8)  +  ik(2W"¥*}  ~WV(?Y  - 

=  -*  k(W*{fV'  +  iW  $<8)')  =  -ikR\*\  (4.54a,  b) 

The  asymptotic  representation  of  as  Y  — ►  ±oo  has  the  form 

*i8)  =  c'£)y^v  +  d(£)y±-v  +  o(r-*+,/), 

where,  in  particular,  —  D ^  is  given  by  the  right-hand  side  of  equation  (4.40)  but  with 

M8)  replacing  Simple  manipulation  gives 

M*1  =  -2i"lt(1-4‘'1)C°a(,r‘/)  /*,  / dt-Jj(T  -  (l)nj(r  -  _  (l) 

“K  JC  JC 

(4.55) 

and  we  find  that 

/si  64tl+l'Jfcl_,'(l  —  2i/)r(i  +  i/W$  r°°  i  _ 

-  D- - ,r(r--f)r(r;\)-  l  ^Ir  -  2<)»,(r  -  (4.56, 

Matching  with  the  outer,  inviscid  solution  requires  that 

D?»  -  D L*>  =  B(6>5  -  i-4'©o«,,  (4.57) 

where  the  subscript  ocf  denotes  the  cubic  part  of  the  total  jump  which  is  due  to  the  outer- 
complementary-function  terms. 


5  The  Evolution  Equations 


The  four  possible  forms  of  the  evolution  equation  for  B(t)  simply  correspond  to  matching 
the  solvability  condition  (3.9b)  with  each  of  the  four  jump  expressions  (4.29),  (4.43),  (4.45)  and 
(4.56)  in  turn;  for  the  sake  of  succinctness,  we  do  not  explicitly  quote  these  evolution  equations. 
Instead,  it  is  more  beneficial  to  the  discussion  here  to  form  one  ‘composite’  evolution  equation 

=  -2 b[%cJ H (») 

+eV7+4M*i+  “  *"4,'C))v  +  eV_7+6l'(&i+  -  i-4"©,},  (5.1) 

where  the  explicit  forms  of  the  nonlinear  quantities  appearing  on  the  right-hand  side  of  this 
equation  have  been  derived,  and  are  quoted,  in  §4;  again  H  denotes  the  Heaviside  function.  As 
the  composite  equation  (5.1)  contains  the  sum  of  the  four  (individual)  nonlinearities,  it  has  the 
advantage  of  being  valid  for  all  values  of  «,  t  and  /i  that  lead  to  an  unsteady  (weakly  nonlinear 
and  weakly  viscous)  critical-layer. 

In  this  section  we  shall  discuss  the  parameter  ranges  of  validity  of  each  of  the  four  possible 
forms  of  the  individual  ‘base’  evolution  equations  for  B(t),  before  presenting  some  numerical 
results  for  the  ‘Ji ’-nonlinearity  case.  This  section  is  concluded  with  a  discussion  of  the  expected 
solution  properties  of  the  other  (base)  evolution  equations.  Further  conclusions  are  drawn  in 
the  next  section. 


2ih 


(i +n 


dB  Jil2k 


5.1  Parameter  ranges  of  validity 


The  range  of  validity/application  of  each  of  the  possible  base  evolution  equations  is  sum¬ 
marised  in  Figures  5.1a,b,  for  v  <  0  and  v  >  0  respectively;  the  governing  balances  determining 
the  dashed  curves  separating  the  three  parts  of  region  III  follow  immediately  by  comparing  the 
sizes  of  the  four  terms  on  the  right-hand  side  of  equation  (5.1).  As  an  example,  we  note  that 
(when  v  <  0)  the  nonlinearity  proportioned  to  J\  is  equal  in  size  or  larger  than  the  quintic 
nonlinearity  if  e2ji_3+4‘'  >  e4/i_7+6l/  i.e.  if  c  <  /i2-".  We  see  that  in  each  case,  depending  on  the 
relative  sizes  of  «,  c  and  p,  three  of  the  four  base  evolution  equations  are  applicable.  Also  plot¬ 
ted  are  the  expected  evolutionary  paths  of  the  disturbance  for  each  case;  this  shall  be  discussed 
in  more  detail  later  but  at  present  it  suffices  to  note  that,  based  on  our  assumptions  concerning 
viscous-spreading  effects  resulting  in  an  unstable  linear  disturbance  mode  approaching  a  later 
neutral  state,  initially  our  disturbance  will  lie  in  the  bottom  right-hand  corner  of  region  Illb 
in  each  case.  This  partly  justifies  our  relegation  of  viscosity  to  lower  order  effects,  a s  well  as 
indicating  that  the  base  evolution  equations  with  the  Jj-cubic  and  the  OCF-cubic  nonlineari¬ 
ties  deserve  the  first  attention.  In  the  next  subsection  we  present  numerical  calculations  for  the 
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former  case  (v  <  0);  the  evolution  equation  with  OOF-cubic  nonlinearity  (t/  >  0)  is  discussed 
in  the  final  subsection. 


5.2  Numerical  results  for  the  Ji-cubic 


Recall  that  if,  and  only  if,  v  <  0  then  a  cubic  nonlinearity  proportional  to  Jj  is  possi¬ 
ble.  In  fact  it  is  the  largest  nonlinearity  in  the  (composite)  evolution  equation  when  /x  k* 
and  t  C  /x2-l/;  thus,  as  mentioned  in  the  previous  section,  it  corresponds  to  the  first  nonlin¬ 
ear  evolutionary  stage  for  an  unstable  linear  disturbance  whose  growth  rate  is  diminished  by 
viscous-spreading  effects.  In  this  subsection  we  present  some  numerical  results  for  this  case. 
To  ease  numerical  calculations  the  jump  expression  (4.29)  is  transformed  into  kernel  form, 

D^-D^  =  V'2t/)  rdss2'AV  fd(r<jl-2''G((T)B(T-s)B{T-as)B(T-(\+a)s) 

41/^1 —  v)  Jo  Jo 


=  2?(6$  -  rn'i'u 


(5.2a) 


where 


r2(|  —  *^)(i  +  <t)  i- 


G(<0  =  [f?^MT{211  +<,)(1  -  il,)F'  (I  -  I + "4 + 1  -  ^  (TTff)) 

-(1  +  O  +  ff’X*  -  1v)Fx  (§  -  \  +  "■ 2  -  2x;  it;  ° 

,»(1  +  *)(5 "  2*0(3  -  2*0  (7  3  1  o- 

+ - 8(1^5 - F'{2-‘,’2  +  ,,'2+,’’3-2'/'0'(uT))} 

+(l+<r)(l-<r)-1"x 

['  ..t-i-'O  -  -  vt)-^  , .  f  1  1  .3  x  !2\] 

l  dt  (l  +  „i)i+-(l-^i)}+-  (2-<,+a(-2'7  OF (-2-‘/'-4-2'4- 2'a  ‘  JJ' 

(5.24) 

Here  F\(a,  5,  c,  d;  x,  y)  is  the  hypergeometric  function  of  two  variables  (see  Erdelyi,  1953;  Churilov 
&  Shukhman,  1988).  The  evolution  equation  for  this  case  can  then  be  written  in  the  form 

fif}  roo  r\  .  _ 

— —  =  - )\B  +  72  /  ds  s2~*v  I  do  a 1~2',G(o)B{t  —  s)B(t  —  os)B(t  —  (1  +  cr)s),  (5.3) 
dr  Jo  Jo 


where 


= - * r, .  — ,  72  =  ^-Fin - r  and  G(o)  =  — sgn (1  +  2v)G. 


8i/r4(i  -  1/)  cos(tti/)/i 


(5.4a  -  c) 
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Let  us  consider  the  range  —2  <  v  <  0;  in  Figures  5.2a,b  — 71/J1  and  —  72/^1  are  plotted 
versus  v  for  this  range.  These  figures  illustrate  that  both  7x  and  72  are  positive  in  this  range 
(recall  that,  since  we  are  considering  unstable  modes,  Jx  <  0).  The  lower  value  v  =  —  2  chosen 
here  for  this  numerical  investigation  has  no  special  significance,  it  was  principally  chosen  so 
that  we  calculate  numerical  results  for  a  complete  period  of  the  trigonometric  functions  present 
in  the  analytical  expressions.  Note  though  that,  for  their  study  of  the  related  neutral  inviscid 
Gortler  modes  in  a  three-dimensional  boundary  layer,  Blackaby  &  Choudhari  (1993)  find  that 
this  range  of  i/-values  is  appropriate.  Moreover,  it  may  also  be  argued  that  for  increasingly 
negative  values  of  u  (corresponding  to  the  Richardson  number  becoming  increasingly  negative) 
the  now  destabilising  effect  of  stratification  means  that  the  growth  rates  of  the  more  unstable 
linear  inviscid  disturbances  are  just  too  large  for  viscous-spreading  effects  to  damp  significantly, 
thus  rendering  the  weakly-nonlinear  theory  inapplicable. 

It  is  convenient  for  numerical  calculations  to  introduce  a  so-called  ‘logarithmic  time’  (see 
Churilov  &  Shukhman,  1988;  Shukhman,  1991) 

T  -  <5-5a> 

having  set 

B(t)  =  B06(r)e^T,  (5.56) 

where  the  constant  B0  is  chosen  such  that  6(f)  — ►  1  as  r  — ►  —00.  This  is  done  in  order  to 
reduce  the  number  of  parameters  in  the  evolution  equation.  Equation  (5.2)  now  has  the  form 

flh  too 

— ;  =  /  d*K{a)  /  dxx2-4,/e-x6(re-:r/(1+"))6(re-ax/(1+<,))6(re-;r),  (5.6a) 

01  Jo  Jo 

with  initial  condition 

b(T  =  0)  =  1;  (5.66) 

where  the  kernel 

K(°)  =  (5-6C) 

In  Figure  5.3  K{a)  ~  K(a\v)  is  plotted,  versus  cr,  for  a  few  representative  v  values.  It  is 
interesting  to  note  that  K(a)  is  always  negative  (0  <  a  <  1);  this  appears  to  be  the  case  for  all 
—2<v<  0;  cf.  the  viscous-jump  kernel  for  the  v  =  0  (J  =  1/4)  case  considered  by  Churilov  &: 
Shukhman  (1988)  which  changes  sign.  Thus  it  is  possible  to  deduce  all  the  qualitative  results 
of  the  solution  properties  of  the  ^volution  equation  (5.6)  from  results  plotted  in  this  figure  (see 
the  discussion  in  §5.2).  However,  we  still  chose  to  solve  the  integro-differential  equation  (5.6) 
numerically  for  completeness,  as  well  as  to  obtain  actual  quantitative  results. 

In  Figure  5.4,  we  present  the  results  of  a  numerical  solution  of  the  evolution  equation  (5.6) 
for  two  representative  v  values.  The  results  show  that  6(T)  oscillates  with  a  fast  rising  amplitude 
as  T  increases.  The  period  of  these  oscillations  is  seen  to  depend  strongly  on  v:  for  v  values 
corresponding  to  the  smaller  x1~Ave~xK(a)  values  (i.e.  u  =  -0.4)  the  time  T  for  6  to  attain  a 
given  large  value  is  longer  than  for  v  values  corresponding  to  the  larger  x2~4t/e~x K(a)  values 
(i.e.  v  —  —1.4).  In  summary,  the  numerical  calculations  indicate  that,  for  all  0  <  v  <  2,  the 
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amplitude  B(t)  will  oscillate,  with  its  magnitude  rising  sharply  after  a  few  oscillations;  our 
results  indicate  that  no  ‘finite-time’  singularity  of  the  evolution  equation  (5.6)  occurs,  instead 
the  oscillations  become  increasingly  wilder. 


5.3  Discussion 


We  begin  this  discussion  by  continuing  to  consider  the  evolution  equation  (5.6)  with  non¬ 
linearity  (solely)  due  to  J\\  later  in  this  subsection  we  consider  the  other  (base)  evolution 
equations.  In  the  last  section  we  reported  that  our  numerical  computations  show  that  K(cr)  is 
everywhere  negative  (0  <  o  <  1)  for  all  —  2  <  v  <  0.  It  was  remarked  that  from  these  facts 
(and  the  magnitude  of  K  =  K(a\  v ))  it  is  possible  to  deduce  all  the  qualitative  results  of  the 
solution  properties  of  the  evolution  equation  (5.6),  without  needing  to  perform  the  numerical 
integration.  To  illustrate  that  this  is  so,  let  us  consider  the  right-hand  side  of  equation  (5.6). 
Initially  b(T)  >  0  and  thus  the  right-hand  side  of  equation  is  clearly  negative  (as  K(e)  <  0); 
thus  this  equation  tells  us  that  b(T)  will  decrease  in  value  until  such  time  when  the  right-hand 
side  of  equation  (5.6)  is  positive.  The  latter  condition  cannot  be  reached  until  b(T)  becomes 
negative.  However  as  the  right-hand  side  of  equation  depends  on  all  previous  b(T)  values,  there 
is  a  delay  until  the  right-hand  side  of  equation  is  actually  positive;  during  this  period  b(T)  has 
becomes  more  negative.  Once  the  right-hand  side  of  equation  is  positive,  b(T )  grows  until  it 
becomes  positive.  Again  there  is  a  delay  until  the  right-hand  side  of  equation  is  negative,  at 
which  time  b(T)  starts  to  decrease  again,  and  so  on.  It  is  clear  that  the  larger  the  typical 
magnitude  of  the  overall  kernel  x2~4ue~x  K(a),  the  more  effect  we  can  expect  the  nonlinearity 
to  have  i.e.  larger  oscillations  at  earlier  T- values.  Note  that  we  would  not  expect  a  singularity 
to  develop  because  of  the  smooth  mechanism  underlying  the  behavior  of  b(T),  as  described 
above.  In  fact,  it  is  possible  to  show  analytically  that  the  solution  of  the  evolution  equation 
(5.6)  can  only  develop  a  singularity  if  /q  K(cr)d<j  >  0  (see  Churilov  fe  Shukhman,  1988). 

Thus,  disturbances  initially  governed  by  the  evolution  equation  (5.6)  would  soon  become 
so  large  in  magnitude  that  their  evolution  would  move  into  its  second  stage  (region  IIIc  of 
Figure  5.1a)  where  the  largest  (leading  order)  term  in  the  composite  evolution  equation  (5.1) 
is  the  quintic  nonlinearity.  We  note  that  this  result  is  not  at  all  surprising  and  could  have  been 
deduced  as  soon  as  the  scales  and  terms  of  the  critical-layer  expansions  (4.4a-d)  were  deduced: 
as  viscosity  does  not  enter  our  analysis  at  leading  order,  the  base  evolution  equations  due  to 
Ji-nonlinearity,  the  OCF-nonlinearity  and  the  quintic-nonlinearity  will  all  lead  to  unbounded 
amplitude  growth,  whether  by  increasing  disturbance  oscillations  or  by  a  singularity  occurring 
(they  do  not  permit  so-called  equilibrium  solutions).  As  viscosity  does  not  enter  at  leading 
order  it  certainly  cannot  damp  out  the  amplitude  growth.  We  note  that  even  with  stronger 
viscosity  effects,  equilibrium  states  may  still  not  necessarily  be  reached  (see  Goldstein  &  Leib, 
1989).  Moreover,  we  note  that  for  the  marginal  instability  case  considered  by  Churilov  & 
Shukhman  (1988),  the  nonlinearity  due  to  viscosity  permitted  singular  solutions  in  which  the 
amplitude  became  unbounded  at  a  finite  T. 
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6  Conclusion 


In  this  paper  we  have  considered  the  nonlinear  development  of  unstable  disturbances  in 
stratified  shear  flows  where  the  Richardson  number  J  is  less  than  one  quarter  in  value.  Although 
such  modes  are  initially  fast  growing,  we  have  assumed  that  viscous-spreading  effects  result 
in  them  evolving  in  a  linear  fashion  until  they  reach  a  state  where  their  amplitudes  are  large 
enough  but  their  growth  rate  have  diminished  significantly  so  that  amplitude  equations  are 
derived  using  weakly  nonlinear  and  unsteady  critical  layer  theories.  We  have  found  that  four 
different  (base)  integro-differential  amplitude  equations  are  possible,  including  one  due  to  a 
novel  mechanism  (the  lJi ’-cubic  nonlinearity).  The  relevant  choice  of  amplitude  equation,  at 
a  particular  instance,  being  dependent  on  the  relative  sizes  of  the  disturbance  amplitude,  the 
growth  rate  of  the  disturbance,  its  wavenumber  and  the  viscosity  of  the  fluid.  This  richness  of 
choice  of  possible  nonlinearities  arises  mathematically  from  the  indicia!  Frobenius  roots  of  the 
governing  linear  inviscid  equation  (the  Taylor-Goldstein  equation)  not,  in  general,  differing  by 
an  integer.  The  initial  nonlinear  evolution  of  a  mode  will  be  governed  by  an  integro-differential 
amplitude  equations  with  a  cubic  nonlinearity  but  the  resulting  significant  increase  in  the  size 
of  the  disturbance’s  amplitude  leads  on  to  the  next  stage  of  the  evolution  process  where  the 
evolution  of  the  mode  is  governed  by  an  integro-differential  amplitude  equations  with  a  quintic 
nonlinearity.  Continued  growth  of  the  disturbance  amplitude  is  expected  during  this  stage, 
resulting  in  the  effects  of  nonlinearity  spreading  to  outside  the  critical  level,  by  which  time  the 
flow  has  become  fully  nonlinear. 

We  finish  by  mentioning  some  further  points  that  may  be  worthy  of  investigation;  mainly 
these  are  related  to  the  relaxation  of  the  assumptions  made  here  in  this  paper.  Obviously, 
the  inclusion  of  viscosity  at  leading  order  in  the  critical  layer  (i.e.  treating  A h  ~  0(1))  would 
answer  the  questions  posed  in  the  previous  section  as  to  whether  it  can  damp  down  the  unlimited 
growth  in  the  magnitude  of  the  disturbances,  allowing  equilibrium  solutions  to  exist.  Another 
obvious  extension  of  this  work  is  to  three-dimensional  disturbances;  then  it  may  be  possible 
to  consider  interactions  between  two  or  more  disturbances.  Comparision  of  our  theory  with 
experiments  or  large-scale  numerical  simulations  is  necessary  at  some  point.  As  mentioned 
in  the  introduction,  as  well  as  geo-physical  flow  applications,  the  theory  has  applications  to 
aerodynamical  flow  situations.  In  fact,  in  a  related  paper  Blackaby,  Dando  &  Hall  (1993) 
apply  the  ideas  contained  in  this  paper  to  the  problem  concerning  the  nonlinear  evolution  of 
inviscid  Gortler  vortices  in  a  three-dimensional  boundary-layer;  their  problem  is  more  complex 
than  that  for  the  model  stratified  flow  considered  here,  as  the  longitudinal  vortex  problem 
is  necessarily  three-dimensional.  The  theory  developed  in  this  paper  can  also  be  applied  to 
a  nonlinear  study  of  the  inviscid  vortex  instabilities  in  the  three-dimensional  boundary-layer 
flow  above  a  heated  plate. 
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Appendix  A 


This  appendix  contains  some  additional  results,  concerning  the  outer  inviscid  problem  for 
the  second  harmonic,  needed  for  the  case  v  >  0.  It  is  merely  a  generalisation  of  the  unpublished 
Appendix  B  of  Churilov  &  Shukhman  (1987)  to  the  case  v  ^  0. 

Outside  the  critical  layer,  the  first  term  of  the  second  harmonic,  ip^,  satisfies 


L2if>£)  =  Q2  =  B2±  sinh3^1  |y| 


2 

cosh4  y 


+ 


3«7p 

2  sinh4  y 


(A.  1) 


We  introduce  the  functions  fa  and  /&  which  satisfy  the  homogeneous  equation  L2rp^  =  0,  such 
that 


fa  =  +  0(e"2h'1)),  fh  =  efc'/P+5M(l  +  0(e~2M)),  as  y  -*  ±oo. 


The  solution  of  (A.l)  is  then  written  in  the  form 

4”  =  BUM  +  CM,)),  (4.2) 

where  C±  are  constants  and  /  is  a  particular  solution  of  equation  (A.l)  that  does  not  contain 
e±fc\A  +3|v|  as  y  — »  ±oo.  As  y  — »  ±0, 

/  =  +  0(  y1))  +  «M*ni  +  °(y))  +  +  0(»));  (-4-3) 

where  a  and  0  are  to  be  determined.  Also 

/.*~a.,l||(|)+‘'(l  +  0(y))  +  M»li','(l  +  0(!')).  ««  V  -*  ±0,  (-4.4) 
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where  the  aajb  and  fia,b  are  fully  determined  by  the  definition  of  fa,b  i.e.  they  are  determined 
solely  from  the  outer,  inviscid  problem.  Thus,  we  can  write 


4"  =  l»rlti‘'(1  +  0(4))  +  B\  (1  +  0(»)), 

(A. 5a) 

as  y  — *  ±0,  where 

<72±  =  a  +  C±aa,  =  /?  +  C±(3a. 

(A. 56,  c) 

It  is  easy  to  show  from  these  relations  that 

ft2+  _  (W*j- )(«A  ^“-)a<|te+/S2.)  _  (aj+/a2.)]- 

(A.5d) 

To  evaluate  a  and  /?,  we  subtract  the  singular  term  in  (A. 3)  by  introducing 

f(y)  =  f  2(3  _  2v)  Smh  ^  ^  C°sh  2  y ’ 

(A. 6a) 

so  that 

/  =  a|y|*+'' +  +  •  •  • ,  as  y-+±0, 

(A. 66) 

and  /(y)  satisfies 

W  =  O2  =  2^3  _  8inh  '+2,/  l»l cosh  2  y  (8  tanh2  2/  (!  +  2j/))  • 

Then 

/  =  [fa{y)fb{z)  ~  fa(z)fb{y)]  , 

(A. 6c) 

(>4.7) 

where  the  Wronskian  W  =  f'afb  —  fafl  =  2j /(<*„$,  —  ctb/3a)  =  —2 ky/k2  +  3.  Considering  y  — ♦  0, 
we  see  that 


a- l  ^(ctbfa  -  <*afb)dz  and  0  =  J  ^y(A/0  -  Pafb)dz.  (A.8) 

Appendix  B 

Again  following  Churilov  &  Shukhman  (1988),  in  accordance  with  (4.34),  we  write 

M+  -  M_  =  (Af+  -  M_)n  +  (Af+  -  M.) 02  +  (Af+  -  M.)3U  (B.  1) 

and  after  extensive  calculations  we  obtain 
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(A/+— A/_)n  = 


f|-7*'ibl-3‘'2-f+l'(l+2i/)r2(|  +  i/)  roo  r 


7r2ra  -  v ) 


f  dt  f  dt4  f  dti  f  dt2 

Jo  Jo  Jc  Jc 


j  dt3B(r  —  t  —  t4  —  ti)B(r  —  t  —  t4  —  t2)B{r  +  t  -+-  tx  +  t2  4-  t4  —  t3)B(r  —  t  —  t4  —  t3)x 
(tiht3)-$-‘'tl+u(t3  -t-h-ti-  2t4)-*~,/(t3  -t-h-  t2)-%~v{t3  -h-  t2)^vx 
^3^l(*l  “  h)2  +  ^2(^2  —  ta)2  +  t\t2(2t3  —  ti  —  <2)j  —  tl+2l/t2(ti  +  t2)  5  "x 

(*i  -  t2)*~u(ti  +  t2  +  t3)(2t3  -t\-  *2).f(-2  -  -  2>  4  ~  2’  if))  X 


H(ti  —  t2)H(t3  —  t  — 1\  —  t2  —  2^4), 


(B.2) 


(M+  —  Af_)o2  = 


2l+*'ii-^A:l-3‘'(l  +  2i/)r(|  +  v )  r°° 


nT*a  ~  p) 


roo  roo  r 

/  dt  dt\  /  dt2 
Jo  Jo  Jc 


J  dt3B(r  —  t~  ti)B(r  -  t  —  t2)B(r  +  t  4-  <1  +  <2  —  h)B(r  -  t  —  t3)(tti)^+‘/ x 


(^2^3)” 2  *(£3  —  2t  —  <1  —  <2)  2  ^2) 2  ^(ti  +  ^2)  2  u{2t  + 1\  +  ^2) x 

(t  +  h  +  *2)f(~  -  -  |i  l  -  ?  -  *2)tf(<3  —  2t  -  «»  -  <2),  (B.3) 


/  ,,  „  x  2-fr'ii-*'kl-»(l  +  2l/)2r(i  +  V)  foo  ^  foo  ^  foo  JA  f  ^ 

(Af+  — Af_)31  = - vT'a-v) - Jo  dtJo  dhJo  dt3Jcdt 2 


,  *T2(|  —  u)  Jo  Jo  “‘Vo  “‘VC““ 

Jc  dt4B(r  -  t  -t3  —  ti)B(r  -  t  -  t3  —  ti)B(r  +  2t  +  tx  +  t2  +  t3  -  t4)~B{T  -  t3  -  t4)  x 

{hh)^vt^~vt ~4*~\U  -3 1-  -tx  -t2-  2t3)-2~l'(t1  +  t2)-'-v{tx  -  trf-'x 
(t4  -3 1-  2t3)^u(t4  -  2t3)-$-u(t4  -3 1-  2 tx  -  2t2  -  2t3){2t4  -6 1-  3 tx  -  3 12  -  4i3)x 

(4  -  l|)F(-i  -  -  J  -  5;  §  -  5;  -  h)H(U  -3  t-U-U-  2i3).  (B.4) 
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Figure  2.1.  Drazin’s  neutral  curve,  J  =  k?(  1  —  k2),  with  the  thick  line  denoting  the  part  of 
the  neutral  curve  with  v  +ve  and  the  thin  line  that  part  with  v  -ve. 
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Figure  5.1a  A  diagram  of  the  various  regimes  of  the  critical  layer  for  u  <  0. 

I:  viscous,  steady  critical  layer;  Landau-Stuart-Watson  equation. 

II:  strongly  nonlinear,  equilibrium  critical  layer;  Benney  &  Bergeron  theory. 

Ilia:  unsteady  critical  layer;  largest  term  in  integro-differential  equation  (IDE)  is  cubic  and 
due  to  viscosity. 

Illb:  unsteady  critical  layer;  largest  term  in  IDE  is  cubic  and  due  to  Jx. 

IIIc:  unsteady  critical  layer;  largest  term  in  IDE  is  quintic. 

The  thick  line  on  the  diagram  indicates  the  expected  evolutionary  path  of  the  disturbance. 
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Figure  5.1b  A  diagram  of  the  various  regimes  of  the  critical  layer  for  v  >  0. 

/:  viscous,  steady  critical  layer;  Landau-Stuart-Watson  equation. 

II:  strongly  nonlinear,  equilibrium  critical  layer;  Benney  &  Bergeron  theory. 

Ilia :  unsteady  critical  layer;  largest  term  in  IDE  is  cubic  and  due  to  viscosity. 

17/6:  unsteady  critical  layer;  largest  term  in  IDE  is  cubic  and  due  to  OCF. 

I  lie:  unsteady  critical  layer;  largest  term  in  IDE  is  quintic. 

The  thick  line  on  the  diagram  indicates  the  expected  evolutionary  path  of  the  disturbance. 
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Solution  of  the  J\  evolution  equation  for  v  =  —0.4. 
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